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Abstract 

Charmonium sum rules for pseudoscalar *~ state ij c (IS) are analyzed within 
perturbative QCD and Operator Product Expansion. The perturbative part of the 
pseudoscalar correlator is considered at o? s order and the contribution of the gluon 
condensate (G 2 ) is taken into account with a s correction. The OPE series includes 
the operators of dimension D = 6, 8 computed both in the instanton and factorization 
model. The method of moments in MS scheme allows to establish acceptable values 
of the charm quark mass and gluon condensate, using the experimental mass of r\ c - 
In result of the analisys the charm quark mass is found to be fh c = 1.26 db 0.02 GeV 
independently of the condensate value. The sensitivity of the results to various ap- 
proximations for the massive 3-loop pseudoscalar correlator is discussed. 

PACS: 13.35.D, 11.55.H, 12.38 

1 Introduction 

The concept of Operator Product Expansion (OPE) was applied to QCD sum rules in [1] to 
parametrize the nonperturbative effects. The operators of increasing dimension, constructed 
from quark and gluon fields, or condensates, constitute the OPE series, which is added to 
the perturbative ones. In case of heavy quark correlators the quark condensates are not 
essential and OPE series start from the dimension D = 4 gluon condensate 

/ Q s /-ia r<a, \ 

for which the authors of [1] have obtained the estimation 0.012 GeV from vector charmonium 
sum rules. In [2] the charmonium sum rules were studied in pseudoscalar channel and it was 
predicted the mass of the lowest r\ c state 3.00 ± 0.03 GeV. This result was in contradiction 
with available to that time experimental information. Later measurements found the mass 
of rj c close to 3.0 GeV, which was considered as a triumph of QCD. 

Since then various sum rules were analyzed in many publications 1 in order to obtain or 
specify the value of the gluon condensate. In the recent paper [3] the vector charmonium sum 
rules were reconsidered with af-corrections of the perturbative series and a s -corrections to 

1 See [3] for the list of publications 
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the condensate contribution, and up to date experimental data. The analysis of [3] resulted 
to the gluon condensate 0.009 ± 0.007 GeV 4 and c-quark mass m c (m c ) = 1.275 ± 0.015 GeV. 
Despite high accuracy of experimental data and c-quark mass determination, the accuracy 
of the gluon condensate remains ~ ±100%, and zero value is not excluded. 

It also seems reasonable to reanalyze the pseudoscalar sum rule, taking into account the 
information obtained in [3]. Now the mass of i] c is known experimentally with high accuracy 
2979.7 ± 1.5 MeV [4], so we invert the problem and find a restriction on the charm quark 
mass and gluon condensate, imposed by this sum rule. A special attention should be paid 
to the correlation between these two values, since a variation of one parameter leads to the 
change of another. 

The sum rule technique goes as follows. The correlator of the pseudoscalar charm currents 
is defined as: 

We define the pseudoscalar current as J p = d^J^, = ey^^c is axial vector current. Within 
the narrow width approximation the imaginary part is: 

1mIF(f + i0) = Tr^T^-m 2 ,) \(0\J P (0)\V)\ 2 (2) 

v 

The sum goes over the pseudoscalar states 2 with J PC = 0~ + . The correlator (1) is quadrat- 
ically divergent, so the dispersion relation requires double subtraction: 

2 . 2 q 4 [°° ImW (s + iO) 



provided the integral in the rhs is convergent, Cq, C\ are unknown constants. In order to 
suppress the contribution of the higher states in (2) as well as continuum contribution, one 
considers the derivatives of the polarization operator in euclidean region Q 2 = — q 2 > 0, the 
so-called moments: 

, 2 , 8tt 2 / d N",^ 8tt f°°lmW(s + i0) 



™ - ^\-dQ 2 ) w ^ = tL l^r^ ds 

8tt 2 v \(o\no)\v)\ 2 U) 

where n > 2. The matrix elements (0 \ J p (0)\i]) are not known experimentally. But if one 
considers the ratio of some two moments at sufficiently high n, the contribution of the lightest 
state Tj c (lS) becomes dominant and 

M n(Q 2 ) _ m 2 Q2 n ^ OQ (5) 
M p n+1 (Q 2 ) ~ m ^ +Q ' n °° (5) 

This property was exploited to predict the mass of r] c in [2, 5]. An essential point was noticed 
in [5]: the QCD corrections to the moments are large at Q 2 = 0, so the sum rules should be 

2 The next to lightest pseudoscalar state r] c {2S) with mass 3654 ± 6(stat) ± 8(syst) MeV was discovered 
recently [6] 
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considered at Q 2 > 0. Moreover, huge contribution of the dimension 8 operators (G 4 ) to the 
moments at Q 2 = [7] becomes tolerable at Q 2 ~ 4m 2 . 

The subject of this paper is a detailed analysis of the sum rule (5). In the next section the 
perturbative and OPE corrections to the correlator (1) are described. Section 2 is devoted 
to the moments both in the pole and MS scheme for the charm quark mass. In the Section 3 
various contribution to the pseudoscalar sum rule (5) are studied in details for typical values 
of the charm mass and gluon condensate. The higher dimension D = 6, 8 gluon operators are 
calculated both in the instanton and factorization model. In the final Section the restriction 
on the c-quark mass and the gluon condensate (aG 2 ) are obtained. 



2 Pseudoscalar correlator in QCD 

In QCD the polarization function (1) consists of perturbative part and operator product 
expansion: 

ip = n PT + n p OPE (6) 

The perturbative part is determined by its imaginary part via dispersion relation (3). The 
imaginary part is parametrized by the coefficient functions R^' p in the expansion by the 
running QCD coupling a(/i 2 ) = a s (fi 2 )/n: 

ImlP PT ( S + *0) = ^ ]Ti^W)aV) (7) 

k>0 

It is simpler to parametrize the functions R^' p in terms of the pole masse m of c-quark. 
The first two terms do not depend on the scale /i 2 . They are known analytically [8]: 



R (0),P 



= -(7-v 2 )+Av (]n^—^--hiv] + 
2 \ 4 3 / 



1-v 2 4, \ 19 + 2u 2 + 3u 4 , 1 + v 8,, 

In- \--(l + v 2 ) 



x 



r , , , l-v\ /3, 1+v , \ , 1 + v 

2Li 2 + Li 2 + -In— \nv In 



1 + v J V 1 + v V2 2 J 1-v 



12 1-v 3 

(8) 



Here and below v = a/1 — TJz, z — s/ (Am 2 ). The function R^' p is usually decomposed into 
the following gauge invariant parts: 

RW'P = C 2 R^' P + C A C F R%\ P + C F Tn l Rf ) ' p + C F TR^' P + C F TRf> p (9) 

where Ca = 3, C F = 4/3, T = 1/2 are group constants and ni = 3 is the number of light 
quarks. The function R^' p comes from the diagram with massless quark loop. It was found 
in [8] and in our normalization takes the form: 



(2) 
P 



(10) 



where the function 8p^ is given by equation (110) in ref [8]. The function R^) ,p comes from 
the diagram with 2 massive quark loops. For s < 16m 2 it contains only the contribution of 
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virtual massive quarks and has the form [8]: 



id 2) ' P = 2vRePj 1 2) - -R^> p In At , s < 16m 2 (11) 
Q 4 m 2 K J 

(2) 

where Pq is second order correction to the pseudoscalar current vertex from the diagram 
with massive quark loop; it is given by equation (169) in ref [8]. For s > 16m 2 the 4-particle 
cut must be included in R^)' p . It is given by the double integral, eq. (97) in ref [8], which 
cannot be taken analytically. Here, however, the total function R^' p can be replaced by its 
high-energy asymptotic, available to the terms m 8 /s 4 in [9]. 

The functions R^' P and R$a correspond to the diagrams with single massive quark 
loop and various gluon exchanges. They are not known analytically, so one has to use some 
approximations. It turns out, that the moments, computed by the dispersion relation (4), 
are sensitive to the choice of these approximations. The accuracy of the moments becomes 
especially important in MS scheme, where there is a sufficient cancellation between large 
terms (see eq. (26) below), so we describe this point in details. 

The first 8 moments M 2 -M 9 at Q 2 = are known analytically [10]. We will require, that 
the approximations for R^' p and R$a must reproduce these moments with high accuracy 
being substituted into the dispersion integral in (4) 3 . As usual, we shall apply the conformal 
mapping and Pade approximation for the relevant parts of the polarization function IP and 
take the imaginary part after then, see Appendix A for details. Although such approxima- 
tions are constructed so that they reproduce low-g 2 expansion of the polarization function, 
they do not give exact values of the first 8 moments at Q 2 = 0, computed by taking the 
dispersion integral in (4). Indeed, the Pade approximations have extra poles away from the 
cut z = [l,oo) and, strictly speaking, the dispersion relation (3) is not valid for them. The 
approximated formulas for R^' p and R^a ■> usec ^ m this paper, are given in the equation 
(46) of Appendix A. 

The last term in (9) R^' p is the so-called singlet part with 2 triangle quark loops. This 
part contains the 2-gluon cut, which is proportional to the 2-photon decay width of the 
pseudoscalar boson. It is known analytically [12, 13]: 



where 



Ri 2), P = TCf \J^2L (12) 



/.^-i|>| 1 -M 1 - I ),i;f jl) ,;;:! 



The contribution of purely gluonic states to the heavy quark current correlators was discussed 
in details in [14] (3-gluon state in case of vector currents). In the narrow width approximation 
(2) only charmed states are taking into account. Since the 2-gluon state is not associated 
with any charmonium state, we subtract R^J' P in the dispersion relation (3) and take the 
integral from s = Am 2 , in accordance with suggestion of [14]. The approximation for R^' p 
without the 2-gluon cut is given in the eq (46) of Appendix A. 

3 We found the approximations proposed in ref [10], eqs (39), (40) having insufficient accuracy to satisfy 
these requirement. 
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The leading in a s order operator series Hq PE (6) for the heavy quark correlator has been 
computed in [15] up to operators of dimension D = 8. This series can be compactly expressed 
in terms of Gauss hypergeometric functions: 



O 



2k-l 

E 



,(p.i) 



fc>2 j 



( 4m 2 )fc -2 ^ fc , (1/2)fc _ i 



2^1 



1, k + i-2 
k - 1/2 



4m 2 



(13) 



where (l/2)^_i is Pochhammer symbol, Of are the operators of dimension D = 2k. For 
the heavy quark correlators there is single operator of dimension D = 4 



O2 = (9 2 G; U G%) 

and 2 operators of dimension D = 6: 



(14) 



(15) 



where gj^ = G a ^ u . u = | X^=« d s Qln^Q- We choose 7 independent operators of dimension 
D = 8 according to [15]: 



of = , of = ((/r 6 ^^) 

' , Of = ((^G^) 



Of 

)(5) 



Of = (g 5 f abc G; u jlj:) , Of = <^ 3 / a6c O^Gt A O c V;Qa > , Of = <£ 4 j';j; ;Qa > (16) 
The coefficients cjff in (13) can be obtained from [15]: 

eg = (0, 1/4, -1/12) , >•'"> 



c 3, i ~~ 



j'V 


1 2 


3 


4 5 


1 


-3/2 


5 


-22/5 1 


2 


1/3 -8/3 


3 


28/15 -4/3 



C 4,i 



j\i 


1 


2 


3 


4 


5 


6 


7 


1 


1/3 


31 


-258 


1580/3 


-365 


70 





2 


1/6 


10 


-111 


946/3 


-4873/14 


144 


-10 


3 


-4/3 


-16 


168 


-596/3 


-106 


140 





4 


1/3 


122 


-946 


5480/3 


-7435/7 


-20 


92 


5 


-10/3 


12 


260 


-2888/3 


7290/7 


-304 


-24 


6 





-22 


152 


-174 


-1108/7 


296 


-96 


7 


2 


-36 


100 


-128 


1314/7 


-208 


72 



(17) 



The a s correction to the D = 4 condensate contribution was obtained analytically in 
[16]. One could differentiate it n times to obtain the moments. However, we prefer to use a 
dispersion-like relation for this correction, constructed in Appendix B, which is convenient 
for numerical calculation of the moments, especially for high n. 
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3 Moments in MS scheme 



At first let us consider the moments in the pole-mass scheme. In QCD the moments (4) 
are expanded by the running QCD coupling a s {ji 2 ). The approximation, used in this paper, 
includes the following ingredients: 1) the perturbative series up to a 2 order, 2) the operator 
series up to dimension D — 8, and 3) a s correction to the (G 2 ) operator contribution. Adding 
all pieces together, we write down the following expression for the pseudoscalar moments: 

2 4 

M£(Q 2 ) = Ys M i k) ' P (Q 2 ) ak ( m2 ) + M n k ' j) ' P (Q 2 ) + 2 Mi° 2 ^(Q 2 )a(m 2 ) , 

k=0 k=2 j 

(18) 

for definiteness the coupling a = a s /n is taken at the scale jj? = m 2 . As discussed in previous 
section, the perturbative moments are taken without 2-gluon cut (12): 



4m2 K s + Q 2 ) n 



The leading order can be expressed in terms of Gauss hypergeometric function: 

M^(Q 2 )= 1 T(n-l) /n-l,n+l _ 
n W) (4m 2 )™- 2 2(l/2) n 2 '\ ra + 1/2 4m 2 ; 



(20) 



where (a) n = r(a+n)/r(a) is Pochhammer symbol. The higher order perturbative moments 
are computed numerically by (19). 

The contribution of the operators to the moments can be easily obtained by differ- 
entiating eq (13): 

2fc-l 



M {Ok,j), P(0 2, = 4 V %\ {k + i-2) n ( n + 1, n + k + % - 2 
3 ^ (4m 2 )"+ fe - 2 fc! (l/2) n+fc _! 2 1 V n + k -1/2 



Am 2 J 



(21) 



The a s -correction to the D = 4 gluon condensate contribution can be obtained by differen- 
tiating eq (53) of Appendix B: 

W j " 3(4m 2 )™ ^ (i - 1)! L (1 + v) n+l Ji (z + v) n+l J 



where y = Q 2 /(4m 2 ), the constants /f and the functions Ff(z) are given in eqs (50) and 
(52) of Appendix B. Notice, that eqs (19)-(22) are applicable for noninteger n also. 

Similarly to the vector case [3], the a s -corrections to the moments are unacceptably large 
in the pole mass scheme and the series (18) is divergent. The pole mass, in fact, is the mass 
of free quark. Since the quarks exist only in form of strongly bounded states, the physical 
meaning of the pole quark mass is rather unclear; it cannot be found from the sum rules 
with a good accuracy. 

Instead of the pole mass one introduces another effective mass parameter, to improve the 
convergence of the perturbative series. Authors of [2, 5] used the mass, renormalized at the 
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euclidean point p 2 = —m 2 . In this paper we shall use the most popular choice for today: the 
gauge invariant mass in the modified minimal subtraction (MS) scheme taken at the scale, 
equal to the mass itself fh = fh(fn 2 ). The pole mass m is perturbatively expressed in terms 
of fh: 



m 2 



= I + Y jKn a n (m 2 ) (23) 



m 

n>l 



The 2-loop factor was found, in particular, in [17] while the 3-loop factor was recently 
calculated in [18]: 

*-§ 

K 2 = 28.6646- 2.0828 n t = 22.4162 

K 3 = 417.039 - 56.0871 n t + 1.3054 n] = 260.526 (24) 

We put rii = 3 in the last column. Now we reexpand the moments (18) by the QCD coupling 

a{fh 2 ): 

M n{Q 2 ) = £m<*>*(QV(™ 2 ) + EE ?'^''^) + 2 M^ p (Q 2 )a(m 2 ) , 

k=0 k=2 j 



where Mn^ ,p {Q 2 ) are expressed in terms of M^' P (Q 2 ) [3]: 
M^' P {Q 2 ) = M^ p 

M^' P (Q 2 ) = M^ p -K 1 (n-d/2)M^ p + K 1 (n+l)Q 2 M^f 
M^iQ 2 ) = M^-K 1 (n-d/2)M^ p + K 1 (n+l)Q 2 Mi l lf 



(25) 



+ (n - d/2) 



9 {n + l-d/2)-K 2 



M (0) ' p 



+ {n + i) [R 2 - K\[n + 1 - d/2)] Q 2 M^f 
+ ^(n + l)(n + 2)Q^l 2 p 
M^° k ' j) ' p (Q 2 ) = Mi° k ' j) ' p 

MjWW^Q 2 ) = M^' p -K l (n + 2-d/2)M^ p + K 1 (n + l)Q 2 M^ G + f ) ' p (26) 

where d = 4 is the dimension of the pseudoscalar function U P (Q 2 ), all Mn^' p in the rhs 
are computed with MS mass fh. The series (25) is much better convergent than (18). The 
numerical values of the ratios M^/M^ and M ( ° 2 ^/M^ for Q 2 /(4m 2 ) =0,1,2 and 
n = 2 — 30 are given in the Table 1 of Appendix C. Notice, that the values of are 
approximate; other approximations for may lead to the moments M^ 2 \ which differ 
from the numbers of the Table 1 within 5 — 10%. 

The expansion (25) goes by a{fh 2 ). If one takes the QCD coupling at some another scale 
/i 2 , the function M^' p changes: 

a(m 2 ) -> a(^ 2 ) , M^ P (Q 2 ) -> M^ P (Q 2 ) + M^ P (Q 2 ) (3 In 4 (27) 



m 2 



so that the series (25) is /^-independent at the order a 2 . 
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4 Pseudoscalar sum rule 



It is convenient to define a dimensionless ratio of the pseudoscalar moments: 

Tn{Q ] ~ 4m 2 M p n+1 (Q 2 ) 4m 2 ' n ^ °° ^ 

Theoretical ratio depends on the quark mass m, QCD coupling a s and condensates. But if 
the dimensionless parameters Q 2 /(4m 2 ), (aG 2 ) /(4m 2 ) 2 etc. are fixed, the l.h.s. of (28) does 
not depend on the quark mass rh (in fact, the QCD coupling depends on the scale, which 
itself may depend on m; but this dependence is weak within the range of error of in). So 
one may use the ratio (28) to find the MS charm quark mass rh for given condensates and 
QCD coupling. 

The QCD coupling constant a s is universal value and can be taken from other experi- 
ments. As input parameter, it is convenient to take a s at the r-lepton mass [4]: 

a s (m 2 T ) = 0.33 ±0.03, m T = 1.777 GeV (29) 

Using this value as the boundary condition in the renormalization group equation, the QCD 
coupling can be evaluated at any scale. As argued in [3], the most natural scale for a s is 

H 2 = Q 2 + rh 2 (30) 

Indeed, in the limit Q 2 3> m 2 we come to natural massless choice a s (Q 2 ), while at Q 2 = 
it becomes a s (m 2 ). Later we shall vary the scale (30) to check the stability of results. 

The MS charm quark mass is determined from vector charmonium sum rules with high 
accuracy. The analysis of the moments at Q 2 = with a 2 corrections leads to rh — 1.304 ± 
0.027 GeV in [19] and 1.23 ±0.09 GeV in [20]. The authors of [19] neglected the condensate 
contribution, while in [20] the value (aG 2 ) = 0.024 ± 0.012 GeV 4 was employed. In fact, the 
gluon condensate weakly affects on the mass value. But for the condensate determination the 
mass accuracy is especially important: a small mass variation leads to significant condensate 
change. As noticed in [3], the perturbative a s and a 2 corrections to the vector moments 
M n (Q 2 ) in MS scheme are strongly suppressed for Q 2 /(4m 2 ) n/5 — 1 and n > 5. The 
analysis of [3] at Q 2 / (Am 2 ) = 1,2 allowed to determine the c-quark mass with high acuracy: 

fh(rh) = 1.275 ±0.015 GeV (31) 

independently of the condensate value. (If the condensate is fixed, the error in (31) can 
be reduced even further.) This result is close to the recent lattice calculation rh = 1.26 ± 
OM(stat) ± 0.12(syst) GeV [21]. 

For the mass (31) one gets the ratio (28) 

O 2 

r n (Q 2 ) = ^ + 1.37 ±0.03 (32) 

in the limit n — > oo. Which values of Q 2 / (4m 2 ) are convenient for the sum rule analysis? The 
choice Q 2 = is not appropriate, since the perturbative corrections to the moments are large 
for almost all n, even in MS scheme. Large Q 2 > 12m 2 are also dangerous: in particular, 
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Figure 1: Ratio r n (Q 2 ) for Q 2 = Am 2 a) and Q 2 = 8m 2 b) versus n. Lower shaded curve 
is purely perturbative, upper shaded curve is computed with condensate (aG 2 ) /(4m 2 ) 2 = 
2 x 10~ 4 . Hatched curves include the (G 3 ) operator (upper curve) and (G 3 ) + (G 4 ) operators 
(lower curve) computed in the instanton model (33). The errorband of each curve corresponds 
to the error of a s in (29). 

when one changes the scale of a s in (27), the effective expansion parameter a(3 In (Q 2 /fn 2 ) 
becomes large > 0.5. In what follows we shall use two choices Q 2 /(4m 2 ) = 1,2. 

The theoretical ratios r„(4m 2 ) and r„(8m 2 ) are plotted versus n in the Fig la) and lb) 
respectively. The lower shaded curve is purely perturbative, i.e. for (aG 2 ) = 0. The central 
line of the shaded area corresponds to the central value of a s (29), the errorband covers 
the error of the coupling a s in (29). One sees, that the agreement with (32) is achieved 
within relatively narrow range of n: n ~ 16 for Q 2 = 4m 2 and n ~ 24 for Q 2 = 8m 2 . 
If we look at the Table 1, the perturbative corrections to the moments in MS scheme, as 
well as a s correction to the condensate contribution, are minimal here. For higher n these 
corrections grow rapidly and the perturbation theory cannot be trusted here. For lower n 
the perturbative corrections are also large, and the leading order of the D = 4 condensate 
contribution crosses at some point, so the behavior of the a s -series is rather unclear here. 
Moreover, unknown contribution of i] c (2S) and higher states to the experimental moments 
could be significant for low n. 

Now we consider nonzero D = 4 condensate. As an illustration, let us fix the ratio 
(aG 2 ) /(4m 2 ) 2 = 2 x 10~ 4 , which corresponds to (aG 2 ) 0.008 GeV 4 , close to the central 
value obtained in [3]. The ratio r n with this condensate is shown by the upper shaded curves 
in Fig 1. The ratio becomes higher for nonzero condensate, which tells in favor of lower mass 
of c-quark. At Q 2 = 4rh 2 the ratio is even higher, than (32) for all n. 

Several models were employed to estimate the higher dimension D = 6, 8 operators. Here 
we consider the dilute instanton gas model [22] and vacuum dominance (or factorization) 
model. 

Instanton model. The vacuum configuration is codsidered as a dilute gas of noniter- 
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acting instantons with effective radis p c and concentration n . The radius p c varies within 
0.3 — lfm in the literature. Here we shall use the estimation p c = 0.5 fm obtained in [23]. 
The instanton concentration n is fitted to the D = 4 gluon condensate (aG 2 ) = 32 n . Then 
one obtains the following expressions for the D = 6, 8 gluon condensates (15), (16): 

(of, . . . , Of ) = (4, 8, 3, 4, 0, 8, 0) ^ 2 (33) 

The ratio r n with the operators 03,4 computed by the instanton model (33) is show by 
hatched curves in Fig 1. Upper hatched curve includes (G 3 ) operator only, lower hatched 
curve includes both (G 3 ) and (G 4 ). 

The (G 3 ) operator contribution is small. But the contribution of the (G 4 ) operators is 
large in the region of interest. Obviously the place, where the lower hatched curve crosses the 
perturbative one (lower shaded), the sum rule (28) with the operators (33) is not applicable. 
Here the (G 4 ) contribution exceeds the leading order (G 2 ), and the OPE series diverges. 

It is a demonstration, that the higher order operators are essentialy overestimated in the 
instanton model [23]. Moreover, their values strongly depend on the instanton size p c , which 
is not strictly fixed. For this reason we finish the analysis in the instanton model. The main 
outcome of this analysis is relatively small contribution of the operator (G 3 ), which will be 
ignored in what follows. 

Factorization hypothesis. In the factorization model the (G 4 ) operators are propor- 
tional to ((G 2 )) 2 . The operators with the light quark current j a in (15), (16) can also 

be estimated by the factorization, but their size is much smaller. The operator Of with 
derivatives was taken as Of m M 2 Of in [15], where M 2 m 0.3 GeV 2 characterizes the 
gluon virtuality in the vacuum. Alternatively, one may express this operator as 

Of = 2<^ 4 r fec G^G^)+2 0f 

Since we neglect the operators with we take Of = 20f here (both estimations agree 
in the order of magnitude for typical condensates). Summarizing, we write down the D = 8 
operators as: 

(^■■■^H^'I^H^ 2 ^ 

The accuracy of the factorization is expected to be ~ ^-/N 2 , N c = 3 is the color number. 
(The l/N 2 ambiguity of the D = 8 quark-gluon condensate factorization was explicitly 
demonstrated in [24].) Another version of the factorization, which employs the heavy quark 
expansion, was proposed in [25]. 

The ratio r n with the operators (34) is shown by the hatched curves in the Fig 2 for 
the D = 4 gluon condensate (G 2 ) /(4m 2 ) 2 = 2 x 10~ 4 . Comparing the Figures 1 and 2 one 
sees, that the contribution of the (G 4 ) operators in the factorization model is smaller than in 
the instanton one. It allows to establish certain stability region, where the ratio r n remains 
almost unchanged. This region is clearly visible for Q 2 = 8m 2 : at n = 20 — 26 the ratio is 
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Figure 2: Ratio r n (Q 2 ) for Q 2 = Am 2 a) and Q 2 = 8rh 2 b) versus n in the factorization 
model. The (G 3 ) condensate is neglected, the contribution of the (G 4 ) operators according 
to (34) is displayed by the hatched curve. Other notations are the same as in Fig 1. 

r n = 3.40 ± 0.01, which corresponds to the c-quark mass m = 1.260 ± 0.005 GeV. This mass 
is computed for the condensate (aG 2 ) = 0.008 GeV 4 . 

In the same way the mass can be computed for other values of the condensate. A 
restriction on the charm quark mass for different condensates is calculated in the next section. 



5 Restrictions on the c-quark mass and D = 4 gluon 
condensate 

As the main result of the pseudoscalar charmonium sum rule (5), we may establish certain 
restrictions for the c-quark mass fh for a given condensate (aG 2 ). 

At first, let us neglect the higher dimension operators (G 3 ) and (G 4 ). The calculation goes 
as follows. For a given Q 2 / (Am 2 ) one should establish the range of n, where the perturbation 
theory as well as operator expansion can be trusted. It is reasonable to require, that the 
perturbative corrections may not exceed 30 — 40% of the leading term. The most dangerous 
is the a^-correction to the gluon condensate contribution M„ ' . Keeping in mind typical 
size of the QCD couping a s /7r ~ 0.1, let us impose the restriction \Mn G '^ /Mn G '°^\ < A. From 
the Table 1 we find the following range of n: 

G 2 O 2 
n = 14 - 19 for = 1 and n = 22 - 30 for = 2. (35) 

Am 2 Am 2 

The perturbative corrections to the moments are also tolerable in this region: the first 
correction \M^ ] /M^ ] \ < 2.5 and the NNLO correction \MP /M^ ] \ < 17. Then, we take 
some value of (aG 2 ) /(Am) 2 and find the maximal and minimal value of the ratio r n (Q 2 ) 
within this range of n. From these numbers we find the minimal and maximal values of the 
charm quark mass fh. 
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Figure 3: Charm quark mass m{m) versus (aG 2 ) obtained from the pseudoscalar sum 
rule. Shaded area displays the acceptable region when the (G 4 ) condensates are neglected. 
Hatched region corresponds to the (G 4 ) condensates, computed by factorization model. 
Dashed lines show the region boundaries for 2 alternative choices of the a s scale. 

The results are shown by the shaded regions in the Fig 3. Fig 3a) and 3b) display the 
restrictions, obtained from the sum rule (28) at Q 2 = 4m 2 and Q 2 = 8rh 2 respectively. Since 
unknown higher order in a s moments are discarded everywhere, the results depend on the 
choice of the scale, at which a s is taken. The dark area shows the acceptable region for the 
scale (30). The dashed and dotted lines display the boundaries of the acceptable region, if 
rh 2 is added to or subtracted from this scale. The scale dependence is weaker at Q 2 = 8m 2 . 

It is clear from the Fig 3, that the pseudoscalar sum rule prefers lower values of the gluon 
condensate. In particular, for the mass rh = 1.275 ± 0.015 GeV [3] one obtains the upper 
condensate limit (aG 2 ) < 0.008 GeV 4 . 

Now let us include the higher dimension operators. As follows from the instanton model 
analisys, the contribution of the D = 6 condensate (G 3 ) is small in the region of interest. But 
the D = 8 operators change the ratio r n essentially. At some n their contribution exceeds the 
leading condensate (G 2 ); at this point the OPE series is divergent. Let us require that the 
contribution of the D = 8 operators to the moments may not exceed 30% of the D = 4 
condensate contibution. This requirement further reduces the region of n (35) depending 
on the condensate size. In the factorization model (G 4 ) ~ ((G 2 )) 2 , and the region becomes 
smaller for higher condensate (G 2 ). 

The hatched regions in the Fig 3 display the inclusion of the (G A ) operators in the factor- 
ization model. For (aG 2 ) > 0.005 GeV 4 the D = 8 operators change the ratio r n drastically. 
They compensate the leading condensate and the ratio r n becomes almost independent of 
the (aG 2 ) condensate in the stability region. From the hatched area in the Fig 3 one gets 
the following limits of the c-quark mass: 

m c (m c ) = 1.26 ±0.02 GeV (36) 
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indenpendently on the condensate value. The mass (36) is in agreement with the result 
m = 1.275 ± 0.015 GeV, obtained from the vector charmonium sum rules in [3]. 

If the (G 4 ) operators are included, it becomes rather difficult to obtain certain restrictions 
on the condensate value. However, for large condensate the stability region is narrow, and 
the results become unreliable. In particular, for {aG 2 } > 0.015 GeV 4 and Q 2 = Am 2 there 
is no region of n, where the OPE series looks convergent. This sets the natural limit of the 
condensate value, at which the pseudoscalar sum rule works. 
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Appendix A: Approximations for R^ ,p 

Let us define the dimensionless coefficient functions for the perturbative correlator (6) as 
follows: 



IP(g 2 ) = 4mV £ Tl^( q 2 )a(m 2 ) , U^(q 2 ) = ^f°° da , (37) 



for definiteness we take the QCD coupling at the scale fj, 2 = m 2 and put the constants 
Co = ci = in the dispersion relation (3). The 3- loop function H^' p is decomposed into 5 
gauge invariant parts in the same way as R^' p (9). 

At first we consider the nonabelian part II^ P . Its expansion near z = q 2 /(4m 2 ) = 
until z 8 is available in [10]. Then, as usual, we reexpand this series in terms of the variable 
u>, which naturally appears in the perturbative calculations: 



uj = (38) 

l + v / T^ 

The expansion of the polarization operator in uj has appropriate analytical properties, namely 
the cut at z — [l,oo). In many cases the Pade approximation was proved to have better 
accuracy, than Tailor series. The best results (see the discussion in Section 2) were obtained 
for the Pade approximation [5/2]: 

U {2) ' p (u) - -^—x 

NA \U) ~ lg7r 2 X 

-7.43220cj + 73.5001cj 2 + 2.69248cu 3 - 13.3868c/ - 0.91032cj 5 - 1.61120cj 6 

x n (39) 

1 - 0.18579cj - 0.63849a; 2 v ; 

The accuracy of the Pade approximated abelian part U^' p is worse because of Coulomb 
behavior ~ — z near the threshold. It turns out, however, that the expansion in uj 
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converges faster, if the multiplier 1/(1 — uj 2 ) is separated out: 



q 

uf' p {uj) = [10.8487 c + 90.9348 c 2 + 111. 510 c 3 + 73.3467c 4 

+113.363 uj 5 + 74.8959 c 6 + 114.366 uj 7 + 77.2105 uj 8 + 0(uj 9 ) ] 



UJ 



= -At — l — [10.8487 uj + 90.9348 c 2 + 100.662c 3 - 17.5881c 4 

lOTT 2 1 — UJ 2 L 

+ 1.85289 c 5 + 1.54919c 6 + 1.00228c 7 + 2.31461 c 8 + 0(c 9 ) ] 

Now we construct the Pade approximation, which well reproduces all asymptotic and first 8 
moments from the dispersion relation: 

10.8487 c + 145.611 uj 2 + 507.376 uj 3 + 57.361 c 4 - 565.406 uj 5 + 94.514 uj 6 
X (1 - uj 2 ) (1 + 5.03984 uj - 4.75471 uj 2 ) ^ ' 

Few more work should be done to construct the approximation for the singlet polarization 
function uf' p . Its expansion near z = until z 8 is available in [11]. The singlet correlator 
contains intermediate massless 2-gluon state, so the cut starts from z — 0, the expansion in 
[11] has the terms ~ In (—z) and the conformal mapping procedure (z — > uj) is not applicable 
here. As discussed in Section 2, in our sum rules we use the polarization operator without 
2-gluon cut (19), so the correspondent part of the polarization operator should be subtracted 
from the result of [1 1] : 

^ = ^C F T(- y) [ WWf., y=<!L (4!) 

where the function f(z) is given in (12). The integral from z = 1 to oo is regular at y = 
and can be expanded by y in Tailor series. But the integral from z — to 1 requires special 
care, since it behaves as ~ lny at y — > 0. In order to obtain the expansion for small y, we 
suggest to use the following series for the function |/(<2)| 2 for < z < 1 [26]: 

, 3^ „ (n — 1)! Z-k 1 . 

aranVi ^S 2 ^Ss (42) 

Then we obtain the following expansion: 

Jo z 2 (z + y) y 2 ^ 



where the constants 



OO , , x , fe — 1 



3 ^ (fc-1)! ^1 

n 2 L> k(k - n - 2) (1/2)* ^ j 2 
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can be computed analytically for any n with the help of recursive relations: 

2(71 + 1)! J 2 7 ^ 

fe=l 



+ 



16(n + 2) (n + 2)(5/2) 

r 2 n 1 "+ 1 



2 



(1/2)* vr 2 



fc 2 (A;-l)! 2 



P(n + 2) A; 3 
fe-i 



, ivfj M _ i y (3/2), 

A; 2Z -Vj + i/2 fc 2Z ^i 2 (i-i 



(44) 



Now one obtains regular at y = Tailor expansion of the polarization operator U^' p without 
the 2-gluon cut, applies the conformal mapping and constructs the Pade approximation: 

-7.86155c; + 6.98952c; 2 + 3.24217c; 3 - 1.68013c; 4 - 0.31547c; 5 - 0.00464c; 6 

x n (45) 

1 - 0.05236 c; -0.32689a; 2 v ; 

Eventually we take the imaginary part of (39), (40), (45) and obtain the correspondent 
coefficient functions R^^: 



d(2),P 1 

Kg ~ 



C F T 



d(2),P _ 
NA — 


8tt 

T 


imnS>) + 


R^ P = 


8tt 

T 


imnj ),p ( w ) , 


P(2),P _ 
ll 99 


8vr 

T 
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16 



E 

m 



where uo 



1 + iy/z - 1 
1 - V* 3 ! 



(46) 



Appendix B: a s - correct ion to the condensate contribu- 
tion 

The a s correction to the D = 4 gluon condensate contribution was found in [16]. Let us 
parametrize it by dimensionless function f^' p (z): 



U p OPE (s) = ...+ 



4tt 2 



(47) 



where dots denote the leading order operator contribution (13). Here we construct a 
dispersion-like relation for this function, convenient for numerical calculation of the mo- 
ments. We will follow the method, used in [3] for the vector current correlator. 
The imaginary part is: 



lmf^ p (z + i0) = 



7T 



9Qz 3 v 5 



In 



zv 



j ^+Pf(z)(l-^2m*; + ^m(4z) 







8) 
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where the polynomials Pf{z) are given in the Table 1 of ref [16]. Taking the contour integral 
in z-plane around the cut z = [1, oo), one could write down the following dispersion relation: 



i 

7T 
+ 



hnfW*(z + iO), ^ n 2 f p 11 e" 3/2 



1+e 



Z-t 



^ d-ty 



15473 35, . 
W + 36 ln(8f) 



1=1 

6' 1 / 2 11 e- 1 ' 2 

T^t + 128 (i -ty 



384 1 - i 



(49) 



where e — ► and 



11 

256 



fp 

J2 



69 



fP — __ 

J3 — 



197 



256 6 2304 

To simplify the calculation, we represent the imaginary part (48) in the following form: 



(50) 



1 



-ImfV>*(z + iO) = F p {z) + F p '(z) + -F. 



p" i 



IT 



(51) 



where the functions F*(z) grow not faster than (z — 1) 1//2 at z — > 1 and have appropriate 
asymptotic at infinity. Our choice is: 



F p (z) = 



1 



96 z 2 v 



49 



207 



56 + —z + 552^ 2 + ( — - 60^ + 904^ 2 - 4216^ 3 + 3312^ 4 



F p (z) 



1 , l-v 
x — In + 

ZV 1 + V 

1 



140 76 



4 

5120 



x 



z 2 + 2208^ 3 21nw + -ln(4^) 



96 z 2 v 
280 



11 645 



197 1439 



24 



24 



1 , l-v 
— m 

ZV 1 + V 



H z 21nw + - In (Az] 



131 197 , l-v 
in 



96 zv 



6 12 zv l+v 



(52) 



Then, the r.h.s. of (49) can be integrated by parts twice, all divergent in e — > terms cancel 
and the dispersion relation can be brought to the following form: 



i=l 



[i-ty 



+ 



i {z-t) 



(53) 



Appendix C: a s -corrections to the moments 

Numerical values of the perturbative corrections to the moments M^' p / Mn^' p and the a s - 
correction to the condensate contribution M^ 2 ^' p / M ( £ )2 ' > ' p are given in the Table 1 for 
Q 2 /(4m 2 ) = 0, 1, 2 and n = 2 - 30. The coefficient functions M% are defined in MS scheme 
according to (25, 26). Remind, that the expansion (25) goes by a(fh 2 ). If one takes the QCD 
coupling at another scale, the function Mn^' p must be changed according to (27). 
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Table 1: a s -corrections to PS-moments in MS scheme 



Q 2 

n 


M 




M 









Am 2 


8ffi 2 





Am 2, 


8m 2 





Am 2 


8m 2 


2 


2.357 


0.825 


-0.176 


-1.132 


-11.89 


-16.24 


1.489 


3.506 


3.617 


3 


3.87 


2.976 


2.064 


12.85 


-0.466 


-8.633 


-1.58 


3.161 


3.731 


4 


3.976 


3.894 


3.179 


19.28 


9.077 


-0.018 


oo 


2.486 


3.611 


5 


3.508 


4.309 


3.831 


20.86 


15.63 


7.025 


3.91 


1.47 


3.32 


6 


2.716 


4.442 


4.226 


19.8 


19.83 


12.51 


0.34 


-0.021 


2.88 


7 


1.71 


4.391 


4.454 


17.56 


22.25 


16.68 


-1.955 


-2.424 


2.291 


8 


0.55 


4.209 


4.562 


15.11 


23.33 


19.75 


-3.979 


-7.447 


1.536 


9 


-0.727 


3.927 


4.581 


13.18 


23.42 


21.93 


-5.923 


-31.87 


0.567 


10 


-2.098 


3.568 


4.529 


12.33 


22.76 


23.37 


-7.847 


31.4 


-0.713 


11 


-3.545 


3.146 


4.42 


13. 


21.58 


24.19 


-9.773 


11.7 


-2.508 


12 


-5.057 


2.671 


4.264 


15.54 


20.04 


24.51 


-11.71 


6.907 


-5.312 


13 


-6.623 


2.152 


4.068 


20.24 


18.29 


24.41 


-13.66 


4.376 


-10.66 


14 


-8.237 


1.594 


3.837 


27.37 


16.44 


23.97 


-15.63 


2.6 


-26.79 


15 


-9.892 


1.003 


3.577 


37.13 


14.6 


23.26 


-17.62 


1.16 


338.4 


16 


-11.58 


0.383 


3.29 


49.72 


12.85 


22.33 


-19.62 


-0.105 


27.93 


17 


-13.31 


-0.264 


2.98 


65.29 


11.27 


21.23 


-21.64 


-1.269 


15.2 


18 


-15.07 


-0.934 


2.649 


84.01 


9.935 


20.01 


-23.67 


-2.373 


10.4 


19 


-16.85 


-1.626 


2.3 


106. 


8.893 


18.7 


-25.72 


-3.439 


7.707 


20 


-18.66 


-2.336 


1.933 


131.4 


8.205 


17.34 


-27.78 


-4.479 


5.863 


21 


-20.49 


-3.064 


1.551 


160.3 


7.917 


15.97 


-29.85 


-5.504 


4.447 


22 


-22.34 


-3.809 


1.155 


192.8 


8.075 


14.6 


-31.93 


-6.519 


3.274 


23 


-24.21 


-4.568 


0.745 


229. 


8.718 


13.27 


-34.03 


-7.527 


2.252 


24 


-26.1 


-5.341 


0.324 


268.9 


9.884 


12. 


-36.14 


-8.531 


1.327 


25 


-28.01 


-6.128 


-0.109 


312.8 


11.61 


10.81 


-38.26 


-9.533 


0.47 


26 


-29.93 


-6.926 


-0.552 


360.5 


13.91 


9.728 


-40.38 


-10.53 


-0.339 


27 


-31.87 


-7.735 


-1.005 


412.3 


16.84 


8.76 


-42.52 


-11.54 


-1.114 


28 


-33.82 


-8.555 


-1.467 


468.1 


20.41 


7.929 


-44.66 


-12.54 


-1.864 


29 


-35.78 


-9.385 


-1.938 


528.1 


24.65 


7.251 


-46.82 


-13.55 


-2.595 


30 


-37.76 


-10.22 


-2.417 


592.2 


29.58 


6.74 


-48.98 


-14.55 


-3.311 
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